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O ' The problem of an impurity electron in a quantum ring (QR) in the presence of a radially directed 

strong external electric field is investigated in detail. Both an analytical and a numerical approach 
to the problem are developed. The analytical investigation focuses on the regime of a strong wire- 
electric field compared to the electric field due to the impurity. An adiabatic and quasiclassical 
approximation is employed. The explicit dependencies of the binding energy of the impurity electron 
on the electric field strength, parameters of the QR and position of the impurity within the QR are 
' obtained. Numerical calculations of the binding energy based on a finite-difference method in two 

' and three dimensions are performed for arbitrary strengths of the electric field. It is shown that the 

binding energy of the impurity electron exhibits a maximum as a function of the radial position of the 
^ ' impurity that can be shifted arbitrarily by applying a corresponding wire-electric field. The maximal 

binding energy monotonically increases with increasing electric field strength. The inversion effect 
of the electric field is found to occur. An increase of the longitudinal displacement of the impurity 
typically leads to a decrease of the binding energy. Results for both low- and high-quantum rings 
are derived and discussed. Suggestions for an experimentally accessible set-up associated with the 
GaAs/GaAlAs QR are provided. 
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I. INTRODUCTION 



During the last decade electronic and optical properties of low-dimensional semiconductor structures have been stud- 
ied extensively both experimentally and theoretically. Along with long-known systems like quantum wells, quantum 
^ , wires, quantum dots, superlattices the novel confined structures called quantum rings (QR) attract much attention. 
The QR can be viewed as a cylindrical quantum dot consisting of an axially symmetric cavity. The unique topology of 
, the QR leads to remarkable quantum phenomena. In the presence of an axially directed magnetic field the persistent 
■ current and the oscillations of the electron energy as a function of the magnetic flux (Aharonov-Bohm effect) were 
^ I ' found to occur^. 

^ To common knowledge impurities and/or excitons modify considerably the electronic, optical and kinetic properties 

of low-dimensional structures such as QRs. Also these properties are strongly affected by external magnetic and 
"j—^ • electric fields. Today, an extensive literature is available which traces the effects of a magnetic field on free carriers, 
' excitons and impurity states in the QR (see for example Lin and Guc^ and Monozon and Schmclcher'^ and references 
^ , therein) . At the same time the influence of an electric field on the electronic properties of a QR has attracted much 
I ' less attention. The energy levels of free electrons and the oscillator strengths of the interband optical transitions as 
' O ' a function of the radii of the QR and strength of the in-plane electric field were investigated in ref0. Barticevic et 
^ al} studied theoretically the effect of the in-plane electric field on the Aharonov-Bohm oscillations and the optical 
^ \ absorption in the QR in the absence of impurities and excitons. Effects of the eccentricity and an in-plane electric 
• • . field on the electronic and optical properties of elliptical QRs have been considered in Ref. 0. Recently the infiuence 
. ^ of the in-planc electric field on the persistent current in the QR coupled to a quantum wire was studied^. In contrast 
, to refs. |4«5.,6_7 the effect of the impurity centre on the electronic states in the QR subject to an axially directed 
^ ■ magnetic and radially directed electric fields was taken into account in Ref. 0- The influence of the radial electric fleld 
5^ " on the electron was assumed to be much weaker than that of the impurity centre, magnetic field and confinement. 
However the infiuence of a strong electric field on the impurity states in the QR is certainly of interest. The reason 
for this is that a strong electric field induces a considerable polarization of the spatial distribution of the carriers^. 
Note that this may be used in order to modulate effectively the intensity of photocurrents and emission of light 
from optoelectronic devices based on QR structures. The analogous effect relating to quantum wells was reported 
by Mendez et al&. Since the problem of the impurity electron in the QR in the presence of a strong electric field 
is not addressed in detail in the literature we investigate this problem here. In case the electric field is parallel to 
the symmetry axis of the QR the ring topology is preserved but there is no significant effect on the radial states of 
the QR. For a strong in-plane electric field the ring topology is broken leading to the disappearance of the unique 
ring properties. It is therefore most advantageous to apply a radially directed electric field created by a wire whose 



2 



position coincides with the symmetry axis of the QR. In this case the electric field is directed radially and a strong 
influence on the radial motion is foreseen and, equally important, the topology of the QR is preserved. 

In the present investigation of an impurity center in a semiconductor quantum ring in the presence of radially 
directed electric field a two-fold approach is pursued. First we will perform studies to obtain analytical (approximate) 
solutions of the stationary Schrodinger equation in certain parameter regimes. This elucidates in particular the 
behaviour and properties of the impurity binding for these regimes. Second, we perform a complementary numerical 
investigation that covers all possible cases (position of impurity, strength of the electric field, radii and height of the 
QR). Finally, exemplary, a comparison of analytical and computational results is provided. 

In detail we proceed as follows. Section II contains a general description of the set-up i.e. the quantum ring, 
the impurity and the electric field configuration. Section III is devoted to the analytical investigation providing the 
method, the results and their discussion. Section IV begins with an outline of our computational method followed 
by a comprehensive discussion of our numerical results. Section V provides a (naturally limited) comparison of the 
analytical and numerical results. Section VI contains the conclusions. 



II. THE QUANTUM RING, IMPURITY AND FIELD CONFIGURATION 

We consider a QR formed by the revolution of a rectangle around the z-axis. The plane of the rectangle contains 
the z-axis. The QR is bounded by potential barriers of infinite height at the planes z = ±d/2 and impenetrable 
cylindrical surfaces at internal p — a and external p — h radii . The chosen model corresponds to a hard-wall 
confinement potential. Alternatively, Chakraborty and Pietilainen^ proposed a parabolic ring confinement potential 
determined by the radius of the ring p and by the effective frequency fJ. This potential has been very effectively used 
in studies of QR3iSiiiii24i^. A comparison of the above mentioned potential models is provided in ref. fill 

The position of the impurity center rg is determined by the cylindrical coordinates a < po < b , ipo = and —d/2 < 
zq < +d/2. Additionally a radially directed electric field is provided by the field of a charged wire with linear effective 
charge density A whose position coincides with the z-axis. Furthermore we take the conduction band to be parabolic, 
non-degenerate and separated from the valence band by a wide energy gap. 

In the effective mass approximation the equation describing the impurity electron possessing the effective mass p 
at a position r(p, (p, z) subject to the axially symmetric and radially directed electric field has the form 

l^^fl^ 1. ^\ -^In^ - 

\ 2p\pdp^dp p^dip^ dz^J 27reoe ''a Aneoe [p^ - 2 pp„ cos ip + pi + {z- zq)^]^' 

= E^ip,^,z) (1) 

where e is the dielectric constant. 

By solving this equation subject to the boundary conditions 

'^'{p, (f, z) = Q for p = a, p = b and z = ± d/2 (2) 

the total energy E and the wave function 5* can be found in principle. 



*(p,</j,z) 



III. ANALYTICAL METHOD AND RESULTS 



A. Adiabatic Approach 



We assume that the effects of the lateral (within the x-y plane) confinement and of the electric field on the (bound) 
electron are taken to be much stronger than the influence of the Coulomb field of the impurity center. Under this 
condition the motion of the electron parallel to the z-axis is adiabatically slower than the motion in the x-y plane and 
the cylindrical variables p, f and z can be adiabatically separated. In the adiabatic approximation the wave function 
\[' can be written in the form 

vi/(p,^,z) = ejv,™(p,</')/'^^"(^) (3) 

where the function 

exp(im(^) 

0Af,m(p, <^) = 7= RN,rn{p) (4) 

\2tt 
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describes the lateral motion of the electron of energy Ej_N^rn determined by the radial confinement and the electric field. 
i?jv,m(p) is the Nth radial state function {N — 1,2,...) possessing the angular quantum number m = 0, ±1, ±2, . . .. 
It vanishes at p = a and p = b (see eq.Q). The function /^'"^{z) describes the longitudinal motion parallel to the 
z-axis and satisfies the equation 



£ £ /^^■™'(^) + VnM /(^■™)(z) = W^(^^™)/(^.™)(z) (5) 



with the boundary conditions 



and with the adiabatic potential 



jWm)(±rf/2)=0 (6) 



T. (,^ _ f dp \RN,rn{p)\'' .7^ 

VN.m[Z) — —- / T-pr (1) 

cos if + pI + {z ~ zq)^] 

The binding energy Eb — E^^"^ — E of the impurity is defined as usual by the difference between the energy of 
the free electron in the QR E^^^ = E±N,m + h^ir^P /2 p d? , I = 1,2,... and the energy of the impurity electron 
E = E^N,m + ly^^'™) which yields 

Eb=^^-W^^^-); ^ = 1,2,... (8) 
where the energy of the longitudinal state VF^^'™) is obtained by solving eq.©. 



B. The Quasiclassical Lateral States of the Confined Electron in the Presence of the Electric Field 

Let us first consider only the electron in the quantum ring i.e. we omit the Coulomb term and the kinetic energy 
in z-direction in cq.(^. Substituting 

RN,m{p) = P"^^^ exp {x/2) UN,ra{x) (9) 

where 

p^P2expx; p2 = aexp{-xi) ; xi = ; 

s 

, 2 _ "ipE^N.m . „ _ 2^eA 

K - -2 , S - [W) 

we obtain the equation for the function un„i{x) 

'^N.mi^) - {jn^ + Spj X exp (2x)j UN,m{x) = (11) 

The boundary conditions |(2Jl become u{xi) — u(x2) = 0, X2 — ln(6/a) + xi. Further we consider the cylindrically 
symmetric radial state only and the corresponding m = label will be dropped in the following. 

Since equation (|lll) does not allow for an exact analytical solution for arbitrary magnitudes of the electric field 
(X s a quasiclassical approach will be used. This method was developed originally for a free electron in an unbound 
medium in the presence of a radial electric field in ref. ^3- In case of a positively charged wire (s > 0) attracting the 
electron to the internal surface of the QR the wave function un (x) is given by the equation 

UNix) - f — j Q-'^Hx) cos[f(x) - -J (12) 
where Q{x) = p2 s^^^(— a^)^^^ exp(x) and where 



F{x) ^ / Q{t) dt (13) 
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The parameters k and p2 are determined by the Bohr-Sommerfeld quantization rule 



Q{t)dt^TT{N + 1/2) ; N ^0,1,2,. 



(14) 



The limits of the integration in ea. (|14|l x = xi and x = correspond to the near and far turning points p = a 
and p = p2- Note that the distance p2 determines the region of the localization of the electron density. Under the 
condition —xi S> 1 the integration in ea. ll4|) can be performed explicitly thereby providing the expressions for the 
energy of the lateral motion E±n<^ and the radius of the iVth radial state p2N 



Ej_N = ^ In 



^ ^ (P2N\ 



2p 



\ a / 



P2N 



4n 



1/2 



as well as the function F{x) ifT^ for a:i < a; < 

F{x) 7r(A^+ 1/2) ( 1 - 



— ) {N + 1/2) 



'X expx 



-1/2 



where <i?(t) is the probability integrali^. Eas. ljlSI) . H16|) are valid under the conditions 



— « 1 ; 

P2N 

For the narrow QR satisfying the conditions 

(6 — a) <C a, 6 ; 
the Bohr-Sommerfeld quantization rule 



> 1 



(15) 



(16) 



(17) 



ib-a)^s 



2tt^ {N 



if a 



< 1 



Q{t) dt = Tr{N + 1) ; N ^0,1,2,. 



leads to the energy of the lateral motion 



E±N — 



h'n^ {N + 1)^ 
2^(6 -a)2 



3A(6- 



(18) 



This result coincides completely with that obtained in ref . |3 where the Schrodinger equation was solved by means of 
perturbation theory. The energy E±pf ()18|l is the size-quantized energy level in the 2D quantum well of width {b — a) 
perturbed by the homogeneous electric field with the effective strength X{A7reoe a)^^ . 

For the negatively charged wire (A < 0) repulsion of the electron towards the external cylindrical surface of the QR 
leads to the Bohr-Sommerfeld rule 

' \Q{t)\dt = 7T{N + 1/2) ■ iV = 0,1,2,... 



which yields for the lateral energy E^m the result 



E 



2p 



\s\ In 



2p 



3^(iV + l/2)|s| 
26 



2/3 



(19) 



This result is valid under the condition 



27r^ (TV + 1/2)^ 



< 1 



It follows from above that the energy E±m (|19|I is the sum of the lowest potential energy of the electron positioned 
at p = 6 in the presence of the radial electric field and a perturbatively acting QR confinement. 

Next we consider the positively charged wire (A > 0) and relatively wide QR for which the lateral energy E±m is 
given by ea. H15|l . In parallel with this we assume that the condition 



P2N < flO 



(20) 



where ag = 47reoe ?i^/J.~^e~^ is the Bohr impurity radius, holds. This means that the effect of the electric field on the 
electron considerably exceeds the influence of the impurity center, so that the energy of the lateral motion i?xjv of 
the impurity electron is determined by the right-hand side of eq. H15|l . 
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C. The Binding Energy of the Impurity Electron 

1. The 2D Impurity States 



Under the condition 



(21) 



the states of the impurity electron have two-dimensional (2D) character. In this case the total energy E can be written 
in the form 



E — E±p^ + ^E 



N 



(22) 



where E±]\[ is given by eg. 1)15(1 and where AEn is determined by the matrix element of the Coulomb term on the 
left-hand side of eq.Q in which we neglect the dependence on the coordinate z. The matrix element is calculated 
with respect to the functions ©, H12|) . I|l()|) with the result 



AE^ = -4Ry (^) 

\P2nJ 



1/2 



1 - 



In (po/g) 
2 In ip2N/a) 



(23) 



where Ry = 11^/2^0^ is the impurity Rydberg constant. Ea. (|23|l is valid under the conditions H17|) . H21|l . (|20|l and 
a ^ po ^ P2N- The energy E (I22|l is obtained by shifting the energy of the lateral motion of the free electron in the 
QR E\ M (tT3|l by the amount AE^ towards lower values of the energy. The binding energy i?t = S^at — E where 
E is given by ea. H22l) becomes Ef, ~ —AEjsi where AEjy is defined by ea. (|23ll . 



2. The 3D Impurity States 



For the relatively high QR satisfying the conditions 



P2N < d, ao 



(24) 



the motion of the electron parallel to the z-axis described by eq.© should be taken into account. The analysis of 
cqs-Q-JTI) is based upon the Hasegawa-Howard method developed originally in ref. 16 and was worked out in further 
detail in ref. 0[ The details of the application of this method to the problem of the impurity in the QR without 
field and in the presence of a strong magnetic field can be found in ref. ,3!. It allows us to restrict ourselves to the 
transcendental equation for the quantum number determining the energy of the ground state of the longitudinal 
motion W. 

a. High QR (d > uq) The equation for the quantum mmrber n < 1 determining the energy W — —Ry/ri^ has 
the form 



2C + il;{l-n) + ^ + 
2n 



Inf^ 

\aon 



Mpo) 



Gi(zo) + G2(zo) 







(25) 



where C is the Euler constant (~ 0.577) and where is the psi-function, the logarithmic derivative of the gamma 
function r(a;). The dependencies of the energy W on the radial and longitudinal positions of the impurity center are 
given by the functions 



where yo — ~ In ( | and 

\P2N, 



A(po) J '^'(%^^) (1/2-yo) +J/0 - (j/o/7r)'^^exp(-2/o) at po < P2N 

. yo at Po > P2N 



^1/2, 



(26) 



Gi,2 = r(-n) 

\aon 

respectively. Ea. (|27|l is valid under the condition 

. 2n 



exp 



_d_/ 2zo 
aon V d 



-—\ 4sinh( — ^1 expi )— 2cosh[ — - 

aon J \aan J [ \aon J \aon 



< 1 



(27) 



(28) 
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Note that eq. (|28|l does not significantly limit the displacement of the impurity zq from the symmetric plane of the 
QR (z = 0). For the QR of height d ~ 2ao subject to the radial electric field providing the relationship p2N — 0.4ao 
and for the impurity being positioned close to the mid plane (zq — 0) and the internal surface (po — a , J/o ^ 1), 
equation (|25l) yields for the quantum number n ~ 0.5. Even though the impurity is shifted by a considerable distance 
Zq = d/4: the term on the left-hand side of ca. (|28|l is about 0.3. In principle ea. l25() can be solved numerically for 
arbitrary values of the height of the QR d > ao, and the impurity position po, d/2 — \zo\ 3> P2N- However, the 
explicit dependencies of the energy W on the above mentioned parameters can be found for the limiting cases of small 
displacements zq from the symmetric plane of the QR (2zo/d ^ 1), for a < p^ <^ p2N and a maximum (po — 6 S> P2n) 
shift Po of the impurity from the symmetric axis p = 0. 

For small displacements Izoj a-o'^ ^ 1 the dependence of the quantum number n(zo) < 1 and the energy W of the 
ground state as a function of the displacement Zq can be found explicitly from eg. 1251) with the result 



Ry 

W{zo,po) = 2" 



l-2r(l-ni) 



2zo \^ f d, 



2ni 



exp 



d 



aoTLi 



(29) 



where ni is the solution to ea. H25(l for zq = and for any radii po. 

The effect of the radial displacement po is described by the function A(po) (|26|) . For p = p2N we have A = and 
with decreasing (increasing) po the function A(po) increases (decreases) towards the internal (external) boundary of 
the QR. We obtain from ea. lf^ 



A(po) = 



2 

P2N- 



PO 



P2N 



for (po — a)^a 
for |po-p27v|<P2Ar 



(30) 



For small radial displacements from po = p2Ar for which 4nA <C 1 (see below ea. (|30|l ') ea.()25|l gives the approximate 
expression for the quantum number n and then for the energy of the ground state W 



W{zo,A) 



Ry 



1 



4^2 
1 + 2^2 



A(po) 



(31) 



where n2 is the solution to eg. 1)251) for A(p2Ar) ~ and any positions zo- 

In the logarithmic approximation {p2N l^o 1; I In {p2N lo-o) I ^ 1) the quantum number n can be calculated from 
eq. 125|l explicitly 



n I \ ao J 



A(po) 



Gi(zo) + G2(zo) 



(32) 



It follows from eq.Q that the binding energy Ef, of the impurity electron in the high QR can be written in the 
form 



Eb = Ry 



1 



(33) 



where the quantum number n can be found from ea. (|25|l to give particularly in the above mentioned logarithmic 
approximation the expression H32|l . It enables us to investigate qualitatively the dependence of the binding energy £'{, 
on the internal radius a and height d of the QR, the strength of the electric field s, and the position of the impurity 
center zo,po- We emphasize that the logarithmic approximation is used only for a qualitative analysis. 

b. Low QR (d < ao) For positive energies W = Ry/s^ in eq.© using the procedure presented in detail in ref. |3 
we obtain the transcendental equation for the quantum number s 



(fi{s) + In 



^ ) - A(yo) 
aos / 



1-i 
2 



Gi(zo) + G2(zo) 



= 



(34) 



In ea. H34f) the following notations are employed 

~f , f(s) r 1 



2i t r(is) 



+ 2C ^l + tP{l + is) 



1 

2is 



1 

2i 



r(is) r(-i.s) 
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Gi,2(zo) = r(s) 



ReTyi,,i/2(Ti,2) 



Tl,2 = 



ImMi54/2(Ti,2) ' '''' iaos V2 
Eq. I|34|l can be solved explicitly for s ^ 1 to give for the energy of the ground state W 



w 



where 




C - A(po 



1 - 



4J 
(P 



(35) 



(36) 



Eas. (|35|) . H3t)|l are valid under the conditions H24|) and d ^ ttoo. The energy W (I35|) is the size-quantized ground 
energy level of the electron in the quantum well of width d shifted towards lower energies by an amount of ISW H3t)|) 
associated with the impurity field. Substituting the energy W (|35|) into the right-hand side of eq.® taken for Z = 1 
we obtain for the binding energy Ei, of the impurity electron in the low QR the result Ei, = — AVF, where ISW is 
given by eq.ljSEJ. 



D. Discussion of the Analytical Results 



1. The 2D Impurity States. 



The binding energy Ei, of the 2D states has the form E^ — —AE^ where the correction to the lateral energy caused 
by the impurity attraction is given by eq. H23() . It is clear from eq. ()23|l that the binding energy decreases with increase 
of the internal radius of the QR a. Also the binding energy decreases with increasing the radial displacement pq of 
the impurity center from the internal surface of the QR in the region a <ti po P2N- Since the strong electric field 
concentrates the electron density close to the internal surface p a the greater the distance pQ — a between this 
surface and the impurity center is the less the impurity attraction i.e. is the less the binding energy. For the impurity 
positioned in the region a ^ po <^ p2N the binding energy increases with increasing electric field strength oc s. The 
increasing electric field shifts the electron density distributed between p c:i a and p ~ p2N towards the impurity center 
positioned at po <C P2N- This leads to an increase in binding energy. 



2. 3D Impurity States, High QR. 



The binding energy Ei, of the 3D states in the high QR (d > oq) is provided by ea. H33|) . The corresponding 
quantum number n can be found in principle from ea. H25|l . In the logarithmic approximation the quantum number n 
is determined by eq. 1)32(1 . Since the contribution of the second term on the right-hand side of eq. ((33|l is exponentially 
small compared to the size-quantized energy (~ 1/d^) (see ea.()29|l') the binding energy decreases with increase of the 
height of the QR d > ag. Clearly from eas. (|33|l and (|32() we see that the binding energy increases with increasing 
the electric field. The greater the electric field is the less the effective radius of the lateral motion peff ~ P2N and the 
greater the depth of the one-dimensional potential V;v,o IZl governing the longitudinal motion between the bottom and 
top of the QR. This leads to an increase with respect to the binding energy Ei,. Expressions and l|29|) show that 
the impurity being positioned at the mid-plane z = produces the greatest binding energy. The shift of the impurity 
center Zo from the mid-plane z = causes a decrease of the binding energy Ei,. Narrowing the QR or increasing the 
electric field strength increases the energy shift A£'{,(zo) associated with the displacement of the impurity from the 
plane = in both cases. 

The radial shift of the impurity center pQ from the some intermediate position pom towards the internal p = a and 
external p = b surfaces produces similar effects as those induced by the displacement zq from the mid-plane z = 0. It 
follows from eqs. (|33|) . H32|) . II26|I . (|3U|) that for QRs possessing an external radius b comparable to the impurity Bohr 
radius oq the binding energy Ei,{po) at po = a approximates that of po — b both being determined by the quantum 

number n ~ 21n^^^ . Since the parameter A(/9o) (|30() and quantum number n^^ (|32l) both increase with the 

displacement of the position po from the external surface po = b towards the internal one p^) — a we expect that the 
binding energy Ei,(pq) reaches a maximum at a certain radial position pom- The shift of the binding energy induced 
by a radial displacement po possesses a maximum for the impurity positioned at the plane zq — and decreases with 
increasing displacement from this plane. 
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It follows from the above that the corrections to the binding energy induced by the displacements from, say, the 
circle po ~ P2N, = to the region a < po < p2N , l^ol > can cancel each other. 



3. 3D Impurity States, Low QR 



It is clear that the binding energy Ei, = — AW (|3t)|) increases with a decrease of the height d < oq of the QR and 
with increasing strength of the electric field strength. For small displacements 2zo /d <C 1 from the mid-plane zq = 
we obtain from H36|) 



-In 



)]} 



(37) 



where 



i?bi(0) = -8i?y(^ 



C - A{po) + In 



rl/2 



P2N 



is the binding energy of the impurity positioned at the mid-plane zq — and any radial distances po- The shift of 
the impurity from the point zq leads to a decrease with respect to the binding energy. With increasing electric field 
strength we observe an increase of the shift of the binding energy caused by the displacement zq. 
The dependence of the binding energy on the radial position po can be derived from eq. (|36|l 



E,iA) = E,2{0) + 8Ry (^) A (po) cos^ (^) 



(38) 



where 



E,2iO)^2Ry{^) 



-4 



C 



cos 



d2 



is the binding energy of the impurity positioned at the point pq — p2N for which A(p2n) = and at any plane zq. 
The binding energy Eb (|38|l decreases if the impurity center moves from the certain position pom towards the radial 
boundaries of the QR. The shift of the binding energy associated with the parameter A(po) reaches a maximum for 
the impurity center positioned at the plane zg = and decreases with increasing displacement from this plane. For a 
low QR the corrections to the binding energy induced by the radial {po) and vertical (zq) displacements can be chosen 
such that they cancel each other, i.e. in this case there is no resulting change of the energy for specific shifts from the 
circle po ~ P2N, zo = to the region a < po < P2N, \zo\ > 0. Thus the dependencies of the impurity binding energy 
Eb on the strength of the electric field, the height of the QR d and the position of the impurity within the QR po, zq 
are qualitatively the same both for high and low QRs. 



IV. NUMERICAL APPROACH 
A. Computational Method 

Our numerical approach to solve is a finite difference method, described in detail in Refs. Il8ll9l2fll for two- 

dimensional systems and in Refs. l2lll2a for three-dimensional systems. We have solved cq.(Q in cylindrical coordinates 
(/9, 0, z) in a region 

a < p < b 

< < TT 

-d/2 <z<d/2 (39) 

respecting the boundary conditions i.e. with = on the boundaries in coordinates p and z and with condition 
= for (j). Our computational procedure consists of the following main steps. The nodes of the spatial 

0—0, 7r 

mesh are chosen in the domain Q and the values of the wave function at the nodes represent solutions of the initial 
differential equation JQ). Since the domain is bounded with respect to all three coordinates we can use uniform 
meshes. The nodes of these meshes have coordinates pi = a + {b — a){i — l/2)/Np, (j>j = 7r(j — l/2)/N^, and 
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Zk = —d/2 + d{k — l/2)/Nz, i = 1, . . . , Np, j = 1, . . . , N^, k = 1, . . . , N^. After replacing the derivatives by their 
finite-difference approximations eq.(l) takes a form of a system of linear equations for "if values at the nodes, and 
approximate values of energy can be found as eigenvalues of the corresponding Hermitian matrix. The final values 
for the energy are provided by using the Richardson extrapolation technique for the corresponding results emerging 
from a series of geometrically similar meshes with different number of nodes, i.e. eigenvalues obtained for meshes 
with Np = KNpo, Ncj, — KNcI,q, Nz — KNzo, where K — 1,2,.... Using this approach we achieve a major increase 
of the numerical precision and, in particular, we obtain together with each numerical value a reliable estimate of its 
precision. Typical numbers of mesh nodes used in the present calculations were of order 40^^, i.e. 40 nodes in each 
direction for the thickest (corresponding to maximal values of K) meshes. An important factor affecting the choice 
of values KNpQ, KN^q, and KN^q is the position of the Coulomb center, which does not coincide with the origin of 
the coordinate system. A geometrical similarity of meshes with different K can be achieved only when this center 
has coordinates {a + {b — a)i/Npo, irj/N^o, —d/2 + dk/Nzo). This circumstance affects the choice of coordinates of 
the Coulomb center and sometimes required calculations on meshes containing more nodes than absolutely needed 
for obtaining a satisfactory numerical precision. Along with solutions of eq.l^ in its general 3D form we have solved 
the corresponding equation in two dimensions employing the coordinates (p, (j)) [z — zq). This allows us obtaining the 
binding energy of the electron in the limit d 0. The numerical solution has no additional specific features compared 
to the 3D one. 



B. Numerical Results and Discussion 

We present here detailed results on the binding energy of the electron for two fixed geometries of the quantum 
ring, referred to in the following as (A) and (B) corresponding to realistic experimental parameters (see below). 
Some additional results on the dependencies of this energy on the geometry of the quantum ring are also presented. 
Parameters of the quantum ring (A) are: GaAs is the ring material, GaAlAs is the barrier material, a = 5 nm, 
b = 20 nm, and d = 15 nm. For the quantum ring (B): a = 10 nm, 6 = 40 nm, and c? = 3 nm with InAs as the 
ring material and GaAs as the barrier material. In order to simplify the comparison of the results for rings made 
of different materials we transform the values of parameters into effective atomic units (e.a.u.). Using parameters 
e = 12.56 and /i = O.OGTmg for GaAs and e = 14.5 and fi — 0.023mo for InAs we obtain a — 0.5, b — 2, and d — 1.5 
for ring (A) and a — 0.3, b — 1.2, and d — 0.09 for (B). The case (A) means a ring with its height being comparable 
with the other dimensions, whereas case (B) is a low ring with d <^ a,b. 

In the following we will particularly consider the binding energy Ey^ of the electron as a function of various pa- 
rameters. The binding energy is the difference between the total energy of the electron obtained in our numerical 
calculations and the energy Eq of the electron in the same quantum ring without the impurity center . The latter 
energy consists of two terms 

Eo = i^iio + ^^±0 (40) 

where E^^q and E±o are the energies of the motion in the z direction and of the lateral motion in the plane (</>, p) 
respectively. An exact analytical expression for the first of them is given in eq.(8) and for the ground state looks for 

2 

effective atomic units as E^^q — ■ The second term does not depend on d and can be easily calculated numerically 
for each set of parameters (a, 6, A) by solving a two-dimensional version of eq.(l) for very large values of po- In Figure^] 
we show i?_Lo(A) for several geometries. They include E^q{X) for QR's (A) and (B) as well as two dependences for a 
small inner radius a — 0.125 e.a.u. and different b. For large positive A these functions are near to linear ones owing 
to the concentration of the electron density in a small vicinity of the inner boundary of the quantum ring. In result 
the curves for a = 0.125 e.a.u. and different b coincide for this range of A. The slope of the curves for positive A is 
determined by a and increases with a decrease of this value. For large negative A the functions E±o{X) also have a 
linear form with slopes depending on both a and b. 

In Figure 121 we present E^ as a function of the position (zo,Po) of the impurity center for quantum ring (A). 
The binding energy shows a significant dependence on po and A for zq ^ d/2 whereas for zq being close to ±| the 
dependency on po is much less pronounced. The dependence of E\, on po and A for zq = is presented in Figure |31 
It should be noted that the effect of the radial electric field on the binding energy is much less compared to its effect 
on the total energy. For example, the maximal value of the difference i?b(A = 10) — i?b(A = 0) (for po = 0.875) is 
0.495 e.a.u., whereas the corresponding difference for Eo (or Ej_o) is 14.81 e.a.u. However, the relative change of i?b 
due to the external field is of the order of 100%. 

For A = the binding energy presented in Figure |31 achieves its maximum (at pom) close to the middle point of the 
radial cross-section of the quantum ring. This is due to the fact that the energy of the ground state of a hydrogen-like 
system increases when the system approach a impenetrable potential wall. In the case of a flat inflnite wall the ground 
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FIG. 1: Energy of the lateral motion a free electron in a quantum ring as a function of the charge A of the central wire. Effective 
atomic units are used. 




FIG. 2: Binding energy Eb as a, function of the position of the impurity center for quantum ring (A) A = 4.3 of the central 
wire. Effective atomic units are used. 




FIG. 3: Binding energy E^, as a function of the radial position of the impurity center po for quantum ring (A) for zo = for 
several different values of the linear charge density A of the central wire. Effective atomic units are used. 
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FIG. 4: Binding energy as a function of the radial position of the impurity center po for quantum ring (A) for A = and 
for different displacements zo = 0, zo = d/2 and the corresponding 2D quantum ring (d = 0). Effective atomic units are used. 



state of a hydrogen atom with the nucleus lying on the boundary of the wall is similar to the state 2po of the free 
hydrogen atom and has energy —0.125 a.u. instead of —0.5 a.u. for the ground state of the free atom. Since Eq does 
not depend on the position of the impurity center approaching the boundaries of the quantum ring equally affects the 
value of the total energy and the value of the binding energy Eh- On the other hand, positive values of E'o for spatially 
confined systems (A — 0) increase the binding energy compared to a free impurity center. As a result the impurity 
electron in the quantum ring is more tightly bound for all the parameter values presented in Figure 13 compared to 
the case of an impurity center in a bulk. However, the total energies of the electron are higher compared to the case 
of a bulk. 

The asymmetry of the curve for A = and in particular a higher binding energy at po ^ a compared to po — ^ ^^e 
due to the curvature of the boundaries of the quantum ring. In the case of the inner po = a boundary its curvature 
provides more space for the motion of the electron compared to a corresponding flat wall. In result the motion of the 
electron is less confined and its energy is lower than for the case of a flat wall and is closer to the values for the case 
where the center is far from the boundaries. The opposite curvature of the outer po — b boundary of the quantum 
ring leads to the opposite effect consisting of a decrease of the binding energy for the impurity center near to this 
boundary. 

For both positive and negative A the maxima of the energy curves E'bCpo) exceed that corresponding to A = 0. This 
effect is fully analogous to the quadratic Stark effect, which leads to the decrease of the ground state energy level and 
to the increase of the binding energy of the hydrogen-like impurity electron in the presence of external electric field. 
The shifts of the positions of the maximums to the left hand side for A > and to the right hand side for A < are 
explained by the increase of the binding energy for smaller distances po in case of an attractive positively charged 
wire or for larger distances in case of a negatively charged wire. 

Figure 3 (curves corresponding to A = ±4.3 ) demonstrates that in the impurity QR the inversion effect of the 
electric field occurs. The binding energy Eb{X) changes as the direction of the electric field A changes (+A ^ — A). 
An analogous effect relating to a quantum well structure was studied in Ref. 12^ In contrast to the quantum well in 
which the inversion shift of the binding energy A_Ef, = Eb(-\-X) — £'f,(— A) vanishes for the impurity centre positioned 
at the mid-plane {zq = 0) , the inversion shift in the QR is absent for a certain {po = 1.15) cylindrical surface. 

In Figure^ we compare the binding energies £'b(po) for the quantum ring (A) for zq = 0, zq = d/2, and for a 2D 
quantum ring with the same values of a and b. We observe a much weaker dependence of the binding energy for 
Zq = d/2 on the radial position of the impurity center compared to the case zq — 0. On the other hand, the curve for 
Zq = is much lower in energy and demonstrates a weaker dependence on po compared to the 2D curve. The large 
difference between 2D and 3D binding energies is due to the large value of d for the quantum ring (A) and strong 
confinement provided by the 2D QR. The predominant part of the curve E^ for the two-dimensional impurity is above 
Ei, = 2. The latter is the binding energy for a two-dimensional impurity without external confinement for the motion 
in p-direction i.e. a bulk impurity. The above-mentioned difference reflects the increase of the binding energy due to 
the confinement in the radial p-direction. 

The opposite case of a low quantum ring is shown in Figure |S1 where the corresponding energy curves analogous to 
those of Figure 21 are presented for the quantum ring (B). For A = the behaviour is qualitatively very similar to the 
one observed in Figure 01 with the exception that both 3D energy curves are much closer to the 2D curve than for the 
quantum ring (A). Figure |S1 contains also an energy curve for a non-zero radial electric field, corresponding to A = 10. 
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FIG. 5: Quantum ring (B): same as Figure|l]and a 3D curve for A = 10, zq = 0. Effective atomic units are used. 
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FIG. 6: Binding energy Eh as a function of the central wire charge A for the different impurity positions po, zo in the quantum 
ring (A). Effective atomic units are used. 



One can see that the radial electric field has a minor effect for the low QR (B) compared to the high QR (A). This 
is due to the quantitatively reduced effect of an external electric field on 2D hydrogen like systems compared to 3D 
ones (see Ref l23l25l and references therein). 

The presence of maxima for the energy functions Eh{po) and the shifts of these maxima with changing radial electric 
field strength determine the form of these functions i?b(A) presented in Figure For the impurity being centered 
in z-direction (zq = 0) we present five curves for different positions po in the radial direction. Let us discuss some 
major properties of these curves. It is straightforward to understand the behaviour of the two curves for po = 1 
and Po = 1.25 corresponding to the impurity center positioned near to the middle of the QR in the radial direction. 
These curves have maxima in the vicinity of A = 0. This means that the presence of a strong radial electric field 
(independent of its direction) decreases the binding energy of the electron to the impurity center for these values of 
Po- This decrease of the binding energy originates from a decrease of the electronic density near the impurity center, 
because the radial electric field either attracts the electronic density to the inner boundary of the ring or repels it 
towards the outer boundary. The curves Eh(X) for the two opposite cases when the impurity center is situated directly 
on the inner po = 0.5 or on the outer po = 2 boundary of the QR possess no maxima. They show an increase of the 
binding energy when the radial electric field attracts the electronic density to the impurity center and a decrease of 
i?b when it repels the electronic density to the opposite boundary. If the boundaries would be fiat and infinite in the z 
direction the limiting values of the binding energy both for pQ = 0.5 and A —> -\-oo and for po = 2 and A — ^ — oo would 
correspond to the ground state of a two-dimensional hydrogen atom — 2 (the mutual action of the electric field 
and the impenetrable potential wall would provide a confinement of the electron in the plane). In the QR the outer 
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boundary provides less space for the motion of the electron in the corresponding limit in comparison with the inner 
boundary. In result we may expect that the limiting value of i?b for pq = 0.5, A —^ +00 should be larger than Ei, for 
po = 2, A — > — 00. This is in agreement with the behavior of the corresponding curves in FigurelHl For the two opposite 
limits po — 0.5, A —00 and po — 2, A ^ +00 ii^b(A) becomes asymptotically flat. This is due to the fact that 
when practically all the electronic density is concentrated at the opposite boundary its small redistribution in strong 
fields does not affect (due to different but large A values) the interaction with the impurity center. The energy of this 
interaction depends on the average distance between the electron and the impurity center. Due to the confinement 
of the electron on the inner or outer surface of the quantum ring the latter problem acquires some similarity with 
the electrostatic problem of a charge near a conducting surface. The potential of the interaction of a charge with its 
image in a convex surface is smaller than the interaction of the charge with a mirror charge in a flat or a concave 
surface (compare with ref. l26l). These reasons explain the result E\^{pQ — 0.5, A —00) > i?b(Po = 2, A ^ +00). 

As an example of i?b(A) for an impurity center located at a small distance from one of the boundaries we show 
in Figure El a curve for pg — 0.75 and zq — 0. In agreement with the picture presented above it has a maximum in 
a region of positive A. The binding energy in the vicinity of this maximum is higher than that for maxima of the 
corresponding curves for po = 1 and po — 1.25 because the radial electric field can concentrate a larger electronic 
density at smaller p. Finally we present also £'b(A) for po = 1 and zq = 0.75, i.e. for the impurity being located on 
the top boundary of the QR. The behavior of this curve is similar to curves for zq — 0, but both the absolute values of 
-Bb and their alterations are smaller. Note that in the QR the radial electric field may provide the impurity electron 
to be more stable while in the bulk material the electric field leads to the ionization of the impurity centre. 

The obtained results allow to estimate the values to be expected in an experiment. It follows from Figure 6 that 
the positive shift of the binding energy Ai?f,(A)of the impurity positioned at Zq = and po = 7.5nm in the QR(A) 
caused by a positive electric field A of a charged wire having the linear electron density rie = 5 • 10^ cm^^ amounts to 
AEh ~ 4.5 meV. The relative change is therefore ^^^''(^^ = 0.28. In this electric field a decrease of the binding energy 
AEb ~ —11.6 meV (about 52 percents) occurs if we shift the impurity centre from the symmetric plane po = 10 nm, 
2:0 = to the boundary plane po — 10 nm, zq — 7.5 nm. When the impurity centre moves in the symmetric plane 
(zo = 0) of the QR(A) (see Figure 3) subject to the positive electric field A — 4.3 determined by the hnear electron 
density rig = 4.3 • 10^ cm~^ from the internal boundary po = 5nm to the position po = 10 nm the binding energy 
Ell increases by an amount AEj, ~ 8.1 meV. This is about 58 percents of the binding energy i?t(po) at po = a. In 
this field the inversion shift of the binding energy AEb{\) = Eh{—X) — Eh{+X) for the case zq = 0, po = 15 nm 
is AEf, ~ 7.5 meV. The inversion shift vanishes at po = 11.5 nm. The estimates of the expected values for the 
InAs/GaAs QR can be made accordingly using the parameters of the InAs material, providing values of the same 
order of magnitude as those for the GaAs QR. Thus the obtained effects induced by the radial electric field in the 
impurity QR are detectable in an experiment. 

Concerning the currently available experimental data to our knowledge most of them are related to the persistent 
current occuring in the QR threaded by the magnetic field (see Ref. [g and references therein). The effect of an 
electric field is studied theoreticallji^*^ for the case of a uniform electric field directed parallel to the plane of the 
QR. One of the reason (concerning both theory and experiment) to choose this configuration is that the electric field 
has been treated as a tool acting on the electron states causing in particular the breaking of the axially symmetric 
potential of the QR and mixing the states with different angular quantum numbers. No additional serious experimental 
refinements specifically associated with the QRs are required. However the radially directed electric field considered 
in our paper is capable on the one hand to conserve the axial symmetry of the QR potential and on the other hand 
to modify strongly the impurity states in the QR. Let us briefiy address the question of the experimental realization 
of the additional electric field. A Si wire covered by the exact-position monitored charge being embedded in the inner 
region of the QR offers to be a source of the radial electric field. Particularly these wires of about 10 nm diameter 
are employed in the silicon-based charge-coupled devices2^. An alternative approach would be as follows. The inner 
region of the QR is doped by the shallow impurity centres. Being activated these centres become charged and this 
region can be treated as the source of the radial electric field for the QR. We believe that our results could stimulate 
experiments that contribute to the physics of QRs and their optical-electronic applications. 

V. COMPARISON OF NUMERICAL AND ANALYTICAL RESULTS 

We present in this section an exemplary comparison of our numerical and analytical results. Figure \7\ shows 
numerically calculated energies for a free electron in three different two dimensional quantum rings. These energies 
are given as functions of the logarithm of the positive charge density on the central wire. For a QR with small 
a — 0.125 and very large b — 12.5 the numerical curve for E{X) presented on a double logarithmic scale is not very 
different from a straight line. The formula (15) provides a good approximation for this curve for not too large values 
of A. The discrepancy at large A is due to small values of p2N {p2N < a) given by eq.(15). A major improvement of 
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FIG. 7: Energy of a free electron in three different two-dimensional quantum rings as a function of the positive charge of the 
central wire A. Numerical results and analytical estimations are shown. Effective atomic units are used. 
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FIG. 8: Binding energy of an electron to an impurity center in a two-dimensional quantum ring as a function of the charge of 
the central wire. Numerical results and analytical estimates are shown. Effective atomic units are used. 



the analytical estimation could be achieved by replacing the formula for p2N by 

/47r\^/' 

P2N=i — ] {N + l/2) + a (41) 

The numerically calculated curve for a = 0.125 and 5 = 2 coincides with that for a = 0.125 and b = 12.5 for large 
values of A. As A — > the value of E converges to a finite limit, determined by the energy of an electron in a finite 
2D quantum ring. This energy increases when the ring becomes narrower as can be seen when comparing the curves 
for a = 0.125 and 6 = 2 and for a = 0.5 and 6 = 2. For large A the latter curve shows also approximately a linear 
behaviour but is lower in energy than the curves for a = 0.125. This shift is properly described by eq.(15), but the 
relatively small value 6 = 2 does not allow to obtain estimations for the energy of the free electron for a — 0.5 and 
6 = 2 by eq.(15). For this set of parameters we present in Figure[7|also the energy given by eq.(18). One can see that 
even for such a broad QR this formula gives quite a reasonable approximation to the energy. For (6 — a) ^ a, 6 this 
formula is in very good agreement with the numerical data. 

In Figure |S| analytical estimates given by eq.(23) for the binding energy of an electron in two-dimensional quantum 
rings in the presence of an impurity are compared with the corresponding numerical results. Eq.(23) includes p2N 
values given by eq.(15) and this circumstance restricts the applicability of it for strong fields. Condition p2N ^ oq, 
or P2N ^ 1 in effective atomic units, restrict the applicability of (23) to weak fields. Within these conditions the 
agreement of the estimates and numerical results is good. This is visible for the curves po = 0.25, which fulfill the 
conditions for the validity of eq.(23) in an improved manner compared to the case po = 0.375. 
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VI. SUMMARY AND CONCLUSION 



We have studied analytically and numerically the problem of an impurity electron in a QR in the presence of a 
radially directed external electric field. The two-fold character of our investigation illuminates the physical behaviour 
and properties of our impurity/quantum ring system in a complementary way. The basis of the analytical approach 
is an adiabatic quasiclassical approximation, while a finite-difference method in two and three dimensions was used to 
perform the numerical calculations. For our analytical studies the external electric field is taken to be much stronger 
than the electric field due to the interaction with the impurity. The dependencies of the binding energy of the impurity 
electron on the strength of the external electric field, the parameters of the QR and the position of the impurity center 
within the QR are derived explicitly. 

We have shown that if the height of the QR increases and/or the impurity center displaces from the mid-circle 
(2:0 7^ 0) for any radius of the QR towards the boundary planes zq = ±d/2 of the QR, the binding energy decreases. 
The binding energy reaches a maximum for the impurity positioned at the mid-plane zq = of the QR. For a fixed 
Zq and without the radial electric field the binding energy has its maximum close to the middle point of the radial 
cross-section of the quantum ring. The radial electric field shifts the position of the maximum towards the center of 
the ring in case of a positive charge of the central wire A and in the opposite direction for a negative charge. This 
results in a relatively complicated dependence -Eb(A), which is very different for different distances of the impurity 
from the center of the QR. The maximum value of the binding energy increases with increasing electric field strength. 
The amplitudes of the mentioned dependencies decrease while shifting the impurity towards to the boundary planes 
of the QR. The inversion effect i.e. the change of the binding energy when the direction of the electric field is changed 
to the opposite one is realized in the impurity QR. Estimates of the binding energies for realistic strengths of the 
external electric field and the parameters of GaAs and InAs quantum rings are provided. 

We have demonstrated that a strong radial electric field and ring confinement provide considerable polarization 
phenomena of the impurity states. Strong dependencies of the binding energy of the impurity electron should lead to 
significant changes of transport processes and optical properties of the QRs. The great sensivity of the impurity QRs 
to the radial electric field is useful for its applications in field-effect transistor structures, electrooptical modulators 
and switching devices. 
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